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Abstract 

Predictions for e + e _ — > J/ipi] c from previous studies are made by taking charmonia as a nonrela- 
tivistic bound state and by using nonrelativistic QCD(NRQCD) approach. The predicted cross-section 
is smaller by an order of magnitude than the experimentally observed. We study the process by tak- 
ing charm quark as a light quark and use light-cone wave-functions to parameterize nonperturbative 
effects related to charmonia. The total cross section of e + e~ — > J/iprj c can be predicted, if these wave- 
functions are known. Motivated by studies of light-cone wave-functions of light hadrons, we make a 
reasonable assumption of the forms of light-cone wave- functions. With these light-cone wave-functions 
we can obtain the cross section which is more closer to the experimentally observed than that from 
NRQCD approach. We also discuss in detail the difference between two approaches. 



At Belle the production of double charmonia at e + e collider with y/s = 10.6GeV has been studied. 
The experimental result is given aspQ: 

a(e + e~ — > J/xfjrj c )Bic(r] c — > 4 charged particles) = 33^g ± 9fb. (1) 

Since a branching ratio is smaller than 1, the above experimental result gives a lower bound for the cross 
section. This experimental result is in conflict with theoretical predictions. Theoretical predictions are 
made by taking a charmonium as a bound state of a cc quark. Employing nonrelativistic wave-functions 
for such a bound state one can predict production rates like the one measured at e + e~ colliders. Starting 
from this, the process were studied in [31I3J[S]. From these studies the cross-section is about 2 ~ 5fb. 
Comparing with Eq.(l) the experimentally measured cross section is about an order of magnitude larger 
than theoretical predictions. 

If one takes charm quarks as heavy quarks, a charmonium system can be thought as a bound state 
consisting of a cc quark mainly, in which the c- and c-quark moves with a small velocity. This fact 
enables us to describe such a system by an expansion in the small velocity. A systematic expansion 
can be achieved by using nonrelativistic QCD(NRQCD) 2|. Within this framework inclusive decays 
and inclusive productions of single quarkonium can be studied consistently and rigorously, where a 
factorization of nonperturbative effects can be completed. But it is expected that theoretical predictions 
for charmonia can have large uncertainties in general. There are two important sources of corrections. 
One is of relativistic correction. Because the velocity of a c quark in a charmonium is not very small, the 
relativistic correction is large. This has been shown in different processes studied in |§1 11U1 ITT] . Another 
source of correction is from high order of a s . If one works with NRQCD factorization for the process 
e + e~ — > r} c J/ijj, the charm quark mass m c can not be neglected and large logarithms like ln(m^/s) will 
appear. Those large logarithms can spoil the perturbative expansion in a s and a resummation is needed. 
However it is not expect that these large uncertainties can result in such a large discrepancy. To explain 
the discrepancy it was suggested that the experimental signals for the final state of J/iprj c may contain 
those of double J/ip® and the initial state interaction can enhance the cross section of J/i/jVc^- Indeed, 
the cross-section becomes large if these suggested effects are taken into account. But, the discrepancy 
still remains large. 

It should be noted that for a process involving two quarkonia there is no rigorous theory based on 
NRQCD in the sense of factorization of nonperturbative effects. Although a proof of the factorization is 
not given yet, but it can be the case that the factorization does not hold. The reason is the following: 
NRQCD is only applicable for a quarkonium in its rest frame. For a process involving only one quarkonium 
one can always find such a rest frame through a Lorentz boost. While for a process involving two or more 
quarkonia one can not find in general a frame in which all quarkonia are in rest. Therefore, a complete 
factorization may not be achieved. 

If the center-mass energy yfs is very large, i.e., yfs 3> m c , one can take c-quark as a light quark. Then 
one can use light-cone wave-functions to describe nonperturbative effects of charmonia and a factorized 
form of the production amplitude in terms of these wave-functions and a perturbative part can be 
obtained. Such an approach for exclusive processes was proposed long time ago 8 . Recently, light-cone 
wave-functions have been employed for charmonia to study their production in B-decav |121 I13| and in 
photoproduction|14j . In comparison with the approach based on NRQCD for the process e + e~ — > J/^rj c , 
where the expansion parameter is the velocity, the approach with light-cone wave-function is with the 
expansion parameters as A/yfs, where A is a soft scale and can be Aqcd, "ic and masses of charmonia. 
In this letter we will use this approach to study the process e + e~ — > J/iprj c . We will work at the leading 
order of the expansion. The correction to our result is at order of A/yfs or (A/yfs) 2 . Taking A to be 
the mass of J ftp, one can expect that our result at y/s = 10 GeV has an uncertainty at level of 30% 
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or smaller. Since an expansion in masses is used, there is only one large scale y/s in the perturbative 
part and our prediction will not contain large logarithms like ln(m 2 /s) if one takes higher orders in a s 
into account. However, such large logarithms like ]n(m 2 /s) and Iji(Aq CD /s) will appear in light-cone 
wave- functions. These large logarithms can be resummed with evolution equations of wave-functions, 
like the Efremov-Radyushkin-Brodsky-Lepage evolution equation|15j. 

In principle the ccg components of charmonia will also contribute, if one take light-cone wave- functions 
at twist-3 into account. Although NRQCD factorization may not directly be applicable to the process 
studied here, but it can be used to study light-cone wave-functions of charmonia as those introduced 
below, in which nonperturbative effects can be factorized into NRQCD matrix elements and light-cone 
wave-functions will be proportional to NRQCD matrix elements|16|. The light-cone wave-functions of 
ccg components will be proportional to NRQCD matrix elements containing gluon fields explicitly. It 
is known that NRQCD matrix elements containing gluon fields is small by NRQCD power counting[2], 
although they can be important in some processes because of some mechanism of enhancement. In our 
approach such an enhancement does not exist, we hence neglect contributions from these ccg components. 
We will only consider contributions from cc components. 

We consider the exclusive process: 

e + (pi) + e"(p2) -» 7*(q) -» J/^(p) + Vc{k), (2) 
where momenta are indicated in the brackets. The amplitude can be written as: 

T = ieu{jpi)^v{p 2 )^e^el{p)p a k P nq 2 l (3) 

where e*(p) is the polarization vector of J/tp and J~(q 2 ) is the form factor defined as 

(J/il>(p)Tic(k)\J»\0) = iQ c ee^el{p) Pa kpHQ 2 ), (4) 

where Q c is the charge fraction of c-quark in unit of e. From the above the produced J/tp is transversally 
polarized. With the form factor the cross section can be calculated as: 

3 

a(e + c- -» J/tP Vc ) = A^a 2 Q 2 }^^- 1 1 _ M^l " f +1 dx (i + ^ (5) 
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where x = cos 9 and 9 is the angle between J/tp and e + . We neglect the small mass difference between 
J/tp and r\ c : 

m h ~ « m Vc . (6) 

From Eq.(4) it is easy to see that the helicity in the process is not conserved. For helicity-conserving 
processes one can use the power-counting rule in [T7j to determine the asymptotic behavior of relevant 
form factors when q 2 — > oo. For a process in which the helicity conservation is violated, one can use the 
generalized power-counting rule in 18 to determine the asymptotic behavior of relevant form factors. In 
our case one can obtain that J-(s) ~ s~ 2 when s — ► oo. 

At the leading order of a s the contribution to the form factor comes from four Feynman diagrams, 
one of them is given in Fig.l. The contribution can be written as 

d A k\ d 4 k2 jj-h n„ u ^ \ f ^4„„-ifei-x / 



(J/tP(p) Vc (k)\J»\0) = J ^^-±H^ kl (k u k 2 ,m c ) J d 4 *e- ifc ^(JMp)|c;(x) Ci (0)|0> 

• J Ae-^<r7 c (£;)|c fc (y)q(0)|0>, (7) 
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Fig.l 

Figure 1: One of the 4 Feynman diagrams for the amplitude. 

where the hard part Hf- kl (k±, k2,m c ) is the amplitude for 

7*(g) -»• Ci(fci) + c fc (fc 2 ) + Cj(p - h) + q(fc - k 2 ) (8) 
and all quarks can be off-shell. The hard part reads: 



(p - fci + A; 2 ) 2 a " (p + fc 2 ) 2 -m2 7 ^ 7 (^-9)2-^ r 



+ (fe-fc 2 + y 2lV j H (fc 2 -g) 2 -m c 27 +7 (fc + A;!) 2 - m 2 7 



(9) 



where r/W are indices for color and spin. 

We will use light-cone wave-functions to describe the two charmonia, in which a collinear expansion is 
used. We take an coordinate system in which J/ip moves in the ^-direction and r] c moves in — z-direction. 
The momenta of these charmonia in the light-cone coordinate system read: 

p" = (p+,|i 0,0) + 1£< 

w = (lp fc "> ' ) = l^" (10) 

where Z M = (1,0,0,0) and n M = (0,1,0,0) are two light-like vectors. The transverse directions to the 
light-like directions are denoted with the subscriber !_. In Eq.(7) the quark pair of c(k\)c(p — k\) is 
transited into the J/ip. In general one expects that the dominant contributions for the integration over 
k\ are with k\± ~ A and k± ~ A 2 /2k±' , where A is the soft scale. Similarly the dominant contributions 
for the integration over k 2 are with k 2 ± ~ A and k\ ~ A 2 /2/c 2 _ . Hence for these integrations one can 
expand the hard part around = (kf, 0, 0, 0) and k% = (0, , 0, 0) and in any soft scale: 

-qj^J (z lP + l,z 2 k-n,0) 

-^p J (ziP + l, z 2 k-n, 0) + m c l-g^f 1 ) (ziP + h z 2 k~n, 0) + • • • . (11) 
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This is equivalent to expand the two matrix elements in Eq.(7) along light cones, i.e., the matrix element 
(7] c (k)\ck(y)ci(0)\0) is expanded around x M = (0, x~,0,0) and (rj c (k)\ck{y)ci(0)\0) is expanded around 
= 0,0,0). However, the expansion is not systematic in the sense that the leading term in 
Eq.(ll) will also lead to some contributions which are at the same order of A as those from higher 
orders. The reason for this is clear: A Dirac field like c{x) can be decomposed along a light cone into a 
"good"- and "bad" component. The "bad" component can be solved with the "good" component with 
equation of motion and will lead to a contribution which is suppressed by A in comparison with that 
from the "good" component. This problem can be solved by expanding the matrix elements according 
to twists of operators. Matrix elements of operators with a given twist are Fourier-transformed light- 
cone wave- functions. For those matrix elements with light hadrons, the expansion in terms of light-cone 
wave- functions have been studied in detail^U EOj. One can use the results in ^31201 to write down the 
expansion for quarkonia. Up to twist-3 the expansion is [T!?l |2*U] : 

{Vc{k)\-c k (y) Cl (m = ^U{[l-n l5 ] lk k- J dz 2 e^ k ~y + ^(z 2 )-^-[ l5 ] lk J dz 2 e^ k ~ <j>f (z 2 ) 

+K75bW"/ dz 2 e^ k 'y + ^(z 2 )} + ■■-, 

<J/V>(p)|«(ste(0)|0) = ^{ifJ/^U^iviP^ J dz 1 e i ^-^(z 1 ) 

+fj/^j / ^% l e*_ Lfl (p) J dz^^-^Z!)} + ■■■, (12) 

where the subscriber _L denote the transverse direction to the light-like directions. It should be noted that 
the space-time coordinate x and y in the matrix elements is not on light-cones, in the side of right-hand 
of the above equations the expansion along light-cones is done. The decay constants are defined as: 

<7 ?c (AO|c(0) 7 ^ 75 c(0)|0) = -if Vc k», 

(7/^)15(0)7^(0)10) = fj^mj/^ip), 2fj^m c = fj^mj^. (13) 

In Eq.(12) the numbers in the bracket [• • •] as subscribers indicate twists. The • • • denote twist-4 terms 
and those twist-3 terms which are proportional to the factor 

^-^^A 1 -^)* ' (i4) 

which vanishes when the quark mass m c goes to infinite. We will neglect contributions proportional to 
this factor. The above wave-functions are normalized, i.e., 

/ o 1 ^wi 2 Wi 3 U [2] ,4 3 Ul 3] }W = i. (15) 

With the expansion in Eq.(12) one can calculate the form factor in terms of these light-cone wave- 
functions. If one only takes twist-2 wave- functions and neglects the quark mass m c , the form factor is zero, 
reflecting the fact that the helicity is not conserved. This also implies that the contribution with twist-2 
wave-functions only is proportional to m c and it is at the same order of those contributions in which 
one of twits-3 wave-functions is involved. We keep the contribution of twist-2 by taking the finite quark 
mass into account. It is straightforward to evaluate the form factor in terms of these wave-functions. We 
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obtain: 



8na s (s) 
9 



m 



2m c iMf \jJ2]/ N 



1 



+m J/ ^ [ l ] (z 1 )^(z 2 ) 
2ml 



1 



1 



+ 



;i - z^z 2 
1 



+ 



1 



;i - z 2 )z\ 
i 



+- 



z\(\-z x ) Z 2 (l-Zi) Zi(l-Z 2 ) 2 z\(l-z 2 ) 

1 1 



+ 



z 2 (i- Zl y z\{\-z 2 ) 



i f3l 

It is interesting to note that the wave- function (p a does not contribute at the considered order, as shown 
by our calculation. The correction to our result in Eq.(16) is suppressed by the power of A/y^ or A 2 /s, 
where A can be the QCD parameter Aqcd, the quark mass m c and masses of quarkonia. 

If we know these wave-functions we can give an numerical result for the form factor and hence the 
cross-section. Unfortunately, these wave-functions are not well known at the energy scale we are interested 
in. If the energy scale is very large, these wave-functions approach to their asymptotic form: 



(16) 



4 2 \z) « $\z)^^(z)^Qz(l-z), 



1. 



(17) 



If we take these asymptotic forms of wave-functions to make predictions, we will have end-point sin- 
gularities. These singularities may be regularized by introducing a momentum cut. We regularize the 
end-point singularities by change the integration range as: 



dz\ I dz 2 



l-e 



dz\ 



l-e 



dz 2 . 



(18) 



with s = m c j\fs. For numerical predictions in this letter we take numerical values of parameters as: 



fljc 



10.6GeU, a s (y/s) = 0.1758, m h ~ 3.0GeV, m c ~ 1.6GeV, 



350MeV, 



2m c 



fj/*, /j/^405MeV. 



(19) 



Taking the asymptotic form of the light-cone wave-functions and these parameters we obtain: 

J/iprjc) ~ 1.31fb. (20) 
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It is interesting to note that light-cone wave-functions as defined in Eq.(12) can be calculated with 
NRQCD factorization, in which nonperturbative effects can be parameterized with NRQCD matrix 
elements [16j. It is easy to obtain the leading order results as: 



4 2 \z) = $\z) = i^{z) = 0(z) = ^ ] (z) = 5(z - i). 



Using this type of wave-functions we obtain for the cross section: 

cr(e + e~ -» J/tpVc) 0.706fb. 



(21) 



(22) 



All predictions in the above two cases are too small in comparison with the experimental result in Eq.(l). 
However, the choice of forms of wave-functions in the two cases is not reasonable, because the asymptotic 
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form is only valid when the energy scale goes to infinity and the NRQCD predictions in Eq.(20) are 
only reliable at the energy scale to be m c with possibly large corrections from higher orders in a s and 
relativistic corrections. Here, we have an energy scale as ~ lOGeV, which is not close to m c and 
far from being infinity In general, predictions are sensible to the form of wave-functions. Light-cone 
wave- functions are nonperturbative objects, which can be only determined with nonperturbative methods 
or extracted from experimental results. The most extensively studied one is the wave-function of n and 
p(e.g., see ^lEHlE^)- Motivated by these studies, we can make some models of wave-functions for 
charmonia. 

A model for the twist-2 light-cone wave- function cp n of it was proposed long time ago in |2j , it takes 
the form as 

M*) = M 1 - z ){ 1 + 2 C [ 5 ( 1 " 2z ) 2 " l ] }' ( 23 ) 

with c = 2/3. A study with QCD sum rule gives c = 0.44 at p = lGeV^Hl- ^ should be noted that the 
shape of <p n with these nonzero values of c is dramatically different than the shape with c = 0, i.e., the 
shape of asymptotic form. Motivated by this observation we assume the twist-2 wave-function for rj c to 
be 

(p [2] (z) = 6z(l-*)|l + 0.44- [5(1 -2z) 2 -l]}. (24) 



The twist-3 wave- functions of tt are also studied in 19. . Through a study of recursion relations of moments 
of [ p' n] (z) and with QCD sum rule the form of (j)p(z) is determined at p = IGeV as: 

4 3 '^ = 1 + (0.39 - 2.5pl)C 1 2 /2 (2z - 1) + (0.117 - 4.914/4 )C\ l2 {2z - 1), (25) 

with p n = (m u + nid) 2 /m 2 . C^(x) denotes Gegenbauer polynomials. It should be noted that the terms 
with P-X represent the part of meson-mass correction. This part can be totally different in the case of r] c . 
Based on this fact we assume the wave-function <p p to be the form: 

4f{z) = 1 + (0.39 - 2.5^JC7 2 1/2 (2z - 1) + (0.117 - 4.914p2jC7 4 1 / 2 (2z - 1), (26) 

where we simply replace p n with p Vc and p Vc takes the form as Am 2 a/M 2 c with a free parameter a. With 
similar technics the wave- functions for p is determined at p = IGeV as 20 : 

ip [ l' p] (z) = 6z(l-z)(l + 0.3(5(2z-l) 2 -l)), 

il) [ l' p \z) = l- 1.6248C 2 1/2 (2z-l)-0.413C 4 1/2 (2z-l), (27) 

We assume the wave- functions of J /if) to be: 

^( Z ) = ^l' p \z\ ^(z) = i)f p \z). (28) 

With these wave- functions in Eq. (23,25,27) we obtain the cross-section for differen values of a: 

u(e + e~ — > J/tpr] c ) ~ 7.37fb, for a = 1, 
cr(e + e~ -> J/ipT] c ) ~ 20.1fb, for a = 1.5, 

cr(e + e- -> J /if; 77c) ~ 31.7fb, for a = 1.75. (29) 

The cross-section increases with increasing a. The predicted cross section with a = 1.75 is much larger 
than that predicted with NRQCD approach and is more comparable with experimental result in Eq.(l). 
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It should be emphasized that the forms of used wave-functions are assumed without any solid arguments, 
although it is motivated with those of ir and p. In this model we neglect the evolution effects of light-cone 
wave- functions. The numbers in the light-cone wave-functions we use are calculated with the QCD sum 
rules for n and p at \i = lGeV. The physics reflected by these numbers are actually nothing to do the 
physics of charmonia. However, it shows the possibility to obtain a large cross-section at the order of 
the experimentally observed with our approach. A detailed study of these light-cone wave-functions of 
charmonia is needed to obtain a reliable prediction. 

To summarize: We have studied the exclusive production of e + e~ — ► J/i/)rf c , in which we have taken 
charm quarks as light quarks and used light-cone wave-functions to parameterize nonperturbative effects 
related to charmonia. In comparison with NRQCD factorization, the factorization of our approach may 
be achieved in a more clean way and the perturbative coefficients will not have corrections with large 
logarithms like ln(y / s/m c ) from higher orders, while in the approach of NRQCD factorization, these large 
logarithms exist and call for resummation. The forms of these light-cone wave-functions are known if 
the energy scale is close to m c or is very large. Unfortunately, these wave-functions at the considered 
energy scale, which is not close to m c and far from being very large, are unknown. With a simple model 
of light-cone wave-functions, we are able to predict the cross-section which is at the same order of that 
measured by Belle. But, this model may not represent the physics of charmonia. A systematic study of 
these light-cone wave-functions is required to have a precise prediction. 
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